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ABSTRACT 

Recent observational surveys have uncovered the existence of super-luminous supernovae 
(SLSNe). While several possible explanations have been put forth, a consensus description 
for SLSNe has yet to be found. In this work we study the light curves of eight SLSNe in the 
context of dual-shock quark novae. We find that progenitor stars in the range of 25-35 Mq 
provide ample energy to power each light curve. An examination into the eff'ects of varying the 
physical properties of a dual-shock quark nova on light curve composition is undertaken. We 
conclude that the wide variety of SLSN light curve morphologies can be explained predom- 
inantly by variations in the length of time between supernova and quark nova. Our analysis 
shows that a singular Ha spectral profile found in three SLSNe can be naturally described 
in the dual-shock quark nova scenario. Predictions of spectral signatures unique to the dual- 
shock quark nova are presented. 

Key words: Dense matter - Radiative transfer - Stars: evolution - (Stars:) supernovae: indi- 
vidual: 2005ap, 2006tf, 2006gy, 2007bi, 2008es, 2008fz, PTF09cnd, PTFlOcwr, 2010gx 



1 INTRODUCTION 

The standard astrophysical explanation for a supernova (SN) is that 
the radiated power is generated by energy deposited in an expand- 
ing ejecta through one of three mechanisms: the SN shock trav- 
els through the stellar envelope (Grassberg et al. 1971), radioac- 
tive decay of heavy elements synthesized in the explosion (Ar- 
nett 1982) or a collision with hydrogen-rich circumstellar material 
(GSM) (Chevalier 1982). In 201 1 astronomers working on the Palo- 
mar Transient Factory announced the emergence of a new class of 
SNe that cannot be explained by any of these means (Quimby et al. 
201 1). As described by Quimby et al. (201 1) this new class of SNe 
is at least ten times brighter than a typical type la SN, displays 
spectra with little to no hydrogen, emits significant UV flux over a 
long period of time and has a late stage luminosity evolution that is 
inconsistent with radioactive decay. 

While this hydrogen-poor class of super-luminous SNe 
(SLSNe) is recent admission, the phenomenon of SLSNe as a 
whole has been an open question since the discovery of SN 2006gy 
(Quimby et al. 2007b). Large scale supernovae surveys such as 
the Palomar Transient Factory (PTF) (Rau et al. 2009; Law et al. 
2009), the ROTSE Supernova Verification Project (RSVP, formerly 
the Texas Supernova Search) (Quimby et al. 2005) and the Catarina 
Real-Time Transient Survey (Drake et al. 2009) have uncovered ap- 
proximately ten other SLSNe, some of which contain hydrogen in 
their spectra (SN 2006gy (Quimby et al. 2007b), SN 2008fz (Drake 
et al. 2010)) while others are hydrogen-poor (SN 2005 ap (Quimby 
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et al. 2007b), SN 2007bi (Gal- Yam et al. 2009)). In order to de- 
scribe these powerful explosions our understanding of the evolution 
of massive stars must change. 

1.1 Proposed SLSN Models 

One method being considered to power the radiated energy of some 
SLSNe is a scaled up version of a GSM interaction. A dense, mas- 
sive (~ 20Mo) GSM envelope must enshroud the progenitor star 
at the time of SN explosion in order for the SN shock to generate 
enough energy to power the SLSN light curve (Smith et al. 2008; 
Chevalier & Irwin 2011; Ginzburg & Balberg 2012; Kiewe et al. 
2012). Building such a GSM envelope requires a mass-loss rate of 
M > O.IM© yr~^ over the final 10-100 years prior to SN explosion 
(Moriya et al. 2012; Ginzburg & Balberg 2012). Two possible ex- 
planations for a mass-loss rate on this order are LBV-like mass ejec- 
tions (Smith et al. 2008; Kiewe et al. 2012) or common envelope 
phase of an interacting binary system (Ghevafier 2012). In com- 
mon envelope description, the SN is powered by the inspiral of a 
compact object on the core of the companion star (Chevalier 2012). 
As well as requiring a massive GSM envelope, the conversion ef- 
ficiency of kinetic energy to radiation must be much higher than 
that of a typical SN (~ 1%) in order to power SLSNe light curves 
(Moriya et al. 2012; Ginzburg & Balberg 2012). Alternatively in- 
creasing the kinetic energy of the SN by approximately ten-fold 
would have a similar effect on the energetics as increasing the con- 
version efficiency. The GSM model cannot explain hydrogen-poor 
SLSN, only type II SLSNe with strong hydrogen lines (Kiewe et al. 
2012). 

An alternative description considered for SLSNe is that the ra- 
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dialed energy is converted from the rotational energy of a magnetar 
(Kasen & Bildsten 2010; Woosley 2010) inside a SN envelope. For 
the magnetar model to power the light curve of a SLSN, large PdV 
losses must be avoided by delaying the conversion of the magne- 
tar's rotational energy into radiation. The delayed injection of en- 
ergy into the SN envelope must be isotropically distributed across 
the inner edge of the SN envelope in order to energize the entire 
envelope and generate enough radiated energy to power a SLSN. 
Whether the magnetar model can power a SLSN with the expected 
jet-like (Bucciantini et al. 2009) energy deposition has yet to be 
studied. Although the isotropic magnetar model can provide a de- 
cent fit to the SLSN light curve, it gives no natural explanation for 
the observed emission lines. 

Pair-instability SNe (PISNe) have as well been proposed as 
the underlying energy mechanism for SLSNe. In this scenario an 
extremely massive star becomes prone to y = 4/3 instability, trig- 
gering a SN explosion. Pan et al. (2012) studied the progenitor stars 
for PISN and found that the mass range required for a star to end its 
Hfe as a PISN is ~ 300- lOOOM© and that no star under lOOOM© in 
the low red-shift Universe would be susceptible to PISN. The rise 
time of the PISN light curve is typically much longer than that ob- 
served in the light curves of SLSNe (Kawabata et al. 2009). In order 
to model the hght curve of SN 2006gy Woosley et al. (2007) sug- 
gested that pair-instability ejections had occurred twice and that the 
second ejection caught and crashed into the first ejection. Woosley 
et al. (2007) stated that in order for this collision to occur the kinetic 
energy of the second ejection must be artificially increased. PISN 
was also suggested by Gal- Yam et al. (2009) as a possible explana- 
tion for SN 2007bi because the SN showed no evidence of a CSM 
interaction. In the supplemental material associated with the work 
presented by Gal- Yam et al. (2009) it was noted that essentially 
all the available nuclear energy of the progenitor star must be con- 
verted to kinetic energy in order for the PISN to explain the light 
curve of SN 2007bi. Table 1 summarizes each proposed model's 
explanation for a variety of SLSNe characteristics. 

Ouyed et al. (2002) first suggested that a collision between 
material ejected through a quark nova (QN) and the preceding SN 
envelope could rebrighten the SN (see section 5.4 of Ouyed et al. 
(2002)). This theory was first applied in the context of SLSNe by 
Leahy & Ouyed (2008). In this work we provide the QN as a pos- 
sible universal explanation for SLSNe. In turn we discuss what 
appear to be observational indications that QNe may lay at the 
heart of SLSNe. Observations of the SLSNe studied in this work 
are introduced in section 2. Section 3 summarizes the explosion 
mechanism of the QN as well as the environment which leads to a 
SLSN. Section 4 examines the physics implemented in describing 
the interaction between ejecta of a SN and a QN. Analysis of the 
eff'ects of changing physical parameters on our model light curve 
is undertaken in section 5. In section 6 we compare observations 
of eight SLSNe (SN 2005ap, SN 2006gy SN 2006tf, SN 2007bi, 
SN 2008es, SN 2008fz, PTF09cnd and PTFlOcwr to QNe of dif- 
ferent physical parameters. A discussion of trends found fitting the 
SLSNe and spectral analysis for some targets is presented in section 
7. Finally our conclusions as well as predicted chemical signatures 
of our model are discussed in section 8. 



2 OBSERVATIONS 

For this analysis we have chosen eight SLSNe (SN 2005ap, SN 
2006gy, SN 2006tf, SN 2007bi, SN 2008es, SN 2008fz, PTF09cnd 
and PTFlOcwr) to study in the context of the quark nova (Ouyed 



et al. 2002). Each SLSN target along with the peak magnitude, class 
and proposed models are summarized in Table 2. 

The first SLSN ever observed was SN 2006gy (Quimby et al. 
2007b) which peaked in absolute R band magnitude at approxi- 
mately -22. The spectrum of SN 2006gy is dominated by a broad 
Hof emission line (Smith et al. 2010) and while SN 2006gy was ex- 
ceptionally bright in visible light the event was surprisingly quiet 
in X-rays (Smith & McCray 2007). 

Spurred by the discovery of SN 2006gy, Quimby et al. (2007b) 
found another extremely bright (-22.7 peak absolute R-band magni- 
tude) SN 2005 ap; which remains the brightest SLSN ever observed. 
The spectrum of SN 2005ap shows broad spectral lines (Ha, C ill, 
N III) and similar to SN 2006gy, SN 2005ap was quiet in X-rays 
(Quimby et al. 2007b). 

SN 2006tf was discovered by Quimby et al. (2007a) who noted 
that the spectrum closely resembles that of SN 2006gy. Observa- 
tions of the light curve of SN 2006tf missed the leading edge of 
the SLSN and thus the actual peak magnitude is unclear, however 
the total radiated energy of the SLSN was at least 7 x 10^^ ergs 
(Smith et al. 2008). Smith et al. (2007) noted that the light curve 
of SN 2006tf is characterized by a very slow luminosity decay rate 
(~ 0.01 Mag. day~^). The spectra of SN 2006tf displays a strong 
Ha emission line that remains broad over time (Smith et al. 2008). 
Hot and BJ3 show an interesting evolution in which the red- side of 
the profile remains constants while the blue-side emission wing be- 
comes more prominent with time (Smith et al. 2008). There exists 
a blue-side absorption feature seen in He i /15876 and O i A1114 
that is of comparable width to the blue-side absorption feature of 
Hof (Smith et al. 2008). 

Gal- Yam et al. (2009) discovered the SLSN SN 2007bi and 
identified it as a type Ic SN, noting that there was no sign of a CSM 

interaction. The light curve of SN 2007bi peaks at 21.3 in the R 

band and displays a slow luminosity decay rate (~ 0.01 Mag. day~^) 
(Gal- Yam et al. 2009). As discussed by Young et al. (2010), the rel- 
atively high host galaxy metalicity is inconsistent with a PISN ex- 
planation for the SLSN. The slowly evolving spectra of SN 2007bi 
shows strong oxygen and iron fines (Young et al. 2010). 

Discovered by Gezari et al. (2009), the light curve of SN 
2008es peaks at ~ -22.2 in the R band and shows a fast luminosity 
decay rate (~ 0.042 Mag. day~^). The total radiated energy in UV 
and visible is in excess of 10^^ ergs, however consistent with other 
SLSN, SN 2008es was quiet in X-rays (Gezari et al. 2009). The 
spectrum of SN 2008es is dominated by broad features that lack the 
narrow and intermediate width line emission typically associated 
with a CSM interaction (Miller et al. 2009). The spectral evolution 
of SN 2008es shows that the broad components of the spectral fines 
become more prominent over time (Miller et al. 2009). 

SN 2008fz was discovered by Drake et al. (2010) who found 
the light curve to peak at ~ -22.3 in the V band and shows a similar 
slow evolution to that of SN 2006gy. The spectrum of SN 2008fz 
displays strong Balmer lines that are initially narrow but become 
broad over time, Drake et al. (2010) noted that the Ha emission 
line of SN 2008fz is similar to that of SN 2006gy for the same 
epoch. 

We chose to study PTF09cnd and PTFlOcwr as representa- 
tive members of the new hydrogen-poor SLSN class (Quimby et al. 
201 1). These SLSNe were selected because they display the two ex- 
tremes of light curve morphology for this class, PTF09cnd has the 
brightest and broadest u-band light curve in the class and PTFlOcwr 
the dimmest and narrowest. The luminosity decay rate of both 
PTF09cnd and PTFlOcwr are inconsistent with radioactive decay 
and the spectrum shows no signs of a CSM interaction (Quimby 
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Table 1. Comparison of proposed SLSNe models. 



Model 



SLSNe Property 



PISN 



CSM 



Magnetar 



dsQN 



Energy Mechanism 
Radiation Mechanism 
Progenitor mass (M©) 
X-rays 

Hydrogen in Spectra 

Cause of long-lasting broad lines 

Late stage luminosity 



7 = 4/3 instabihty^ 
Collision of ejecta^ 
313^ - 1250^ 
suppressed^ 
unlikely"^ 
not expected^ 
radioactivity^ 



Binary merger" 
CSM/SN interaction^ 
100+g 
suppressed^ 
necessary^ 
velocity of SN§ 



Rotational energy^ 
Synchrotron*^ 

8*^ - 25' 
not discussed 
not discussed 
not expected'^ 



opaque outer region of CSMP inner bubble'^ 



QN explosion ^ 
QN/SN interaction^ 
20' -40' 
suppressedj 
likelyj 

temperature of inner shelP 
inner shell emission^ 



^Woosley et al. (2007), ''Chevalier (2012), ^Kasen & Bildsten (2010), ^Leahy & Ouyed (2008), ^Chevalier (1982), ^Pan et al. (2012), ^Smith et al. (2010), 
^Weidemann (1977) ^Fryer (1999), JQuyed et al. (2012), ^Blinnikov (2008), ^Chevalier & Irwin (2012), "^Yungelson et al. (2008), ^Quimby et al. (201 1), 

°Kasen et al. (2011), PChevalier & Irwin (2011), 



et al. 2011). PTFlOcwr was also observed in the B band by Pas- 
torello et al. (2010), in which it peaks at ~ -21.2, and alternatively 
is referred to as SN 2010gx. As mentioned in the supplemental ma- 
terial associated with Quimby et al. (2011), the SLSNe PTF09cnd 
and PTFlOcwr are quiet in X-rays relative to the amount of energy 
radiated in the visual bands. 

3 QUARK NOVA 

3.1 Explosive Mechanism 

The existence of up-down-strange (UDS) quark matter as a state 
of matter more stable than hadronic has been theorized for several 
decades (Itoh 1970; Bodmer 1971; Witten 1984). However, Ouyed 
et al. (2002) was first to propose that if a phase transition to UDS 
matter were to occur at the core of a neutron star the explosive re- 
sult would be capable of ejecting the outer layer of the neutron star 
leaving behind a quark star. This violent new pathway for massive 
star evolution is called a QN. It has been shown for massive neu- 
tron stars that either the accretion of SN fall-back material or the 
spin-down evolution can cause the core density to exceed that of 
quark deconfinement, triggering a QN (Staff" et al. 2006). Detailed 
study of hadronic to UDS burning suggested that turbulence in the 
conversion front enhances the burning speed to sonic or supersonic 
velocities which in turn would make the transition from neutron star 
to quark star detonative (Niebergal et al. 2010; Ouyed et al. 201 lb). 
Further investigation of neutrino and photon emission from a QN 
found that 10"^ - 10"^ M© could be ejected with up to 10^^ erg of 
kinetic energy, which means ejecta with a Lorentz factor of 10-100 
(Vogt et al. 2004; Ouyed et al. 2009). 



3.2 Timing 

The length of time between SN and QN explosion (fdeiay) plays a 
crucial role in determining the subsequent evolution of the star. 
When the delay is short (< 8 days) the SN envelope is still dense 
and the energy of the impacting QN ejecta is used up in spallation 
of the inner region of the SN envelope (Ouyed et al. 2011a). This 
leads to the destruction of ^^Ni and the formation of "^"^Ti and re- 
sults in a subluminous SN. The reduced luminosity is due to the 
lack of radioactive decay of ^^Ni (Ouyed et al. 2011a). When the 
time delay is long (on the order of months or longer) the SN enve- 
lope will have become too diff'use to significantly interact with the 
QN ejecta. For massive QN ejecta with long time delays fall-back 
of the ejected material can occur which has implications for y-ray 



bursts, soft y-ray repeaters and anomolous X-ray pulsars (Ouyed 
et al. 2007a,b; Koning et al. 2012). 

This work will focus on the scenario in which a SN is followed 
on the order of weeks by a QN. In this case, referred to as a dual- 
shock QN (dsQN), the expanded SN envelope is bombarded by 
the QN ejecta reheating the SN envelope. The extended size and 
high temperature of the re- shocked SN envelope yields a brilliant 
radiance capable of reproducing the observed luminous peaks of 
SLSNe. As the radiation from the extended envelope fades, an inner 
region of mixed QN and SN material is revealed which can explain 
both the luminosity decay rate of SLSNe as well as curious spectral 
features. 



4 DUAL-SHOCK QUARK NOVA 

The evolution of a dsQN can be considered as three distinct phases; 
delay, shock and cooling. During the delay phase the SN envelope 
expands homologously (v oc r) while the neutron star evolves to- 
wards a QN explosion. For fiducial values the radius of the SN 
envelope will reach ~ 10^^ cm during this phase. The end of the 
delay phase is marked by the detonation of the QN. 

At the beginning of the shock phase the ultra-relativistic QN 
ejecta (Keranen et al. 2005; Ouyed & Leahy 2009) quickly catches 
and slams into the inner edge of the SN envelope. This collision 
creates a shock front that propels through the SN envelope reheat- 
ing it to a temperature ~ 10^ K (Leahy & Ouyed 2008). As the 
shock progresses the inner region of the re-shocked SN envelope 
mixes with the impacting QN ejecta to create a thin shell interior to 
the envelope (referred to as the hot plate in paper 1). 

The end of the shock phase and beginning of the cooling phase 
is defined to be the moment when the shock breaks out of the SN 
envelope. By this time the inner shell will be fully formed and 
slowly coasting inside of the shocked envelope. During the cooling 
phase both the envelope and the inner shell will cool via adiabatic 
expansion and radiative emission. Following the same methodol- 
ogy used in Ouyed et al. (2012) (furthermore to be referred to as 
paper I), for the work presented here we model the light curves of 
dsQNe during their cooling phase. 

4.1 The Inner Shell 

As the QN ejecta ploughs through the SN envelope material builds 
up to form the inner shell. The mass (M^h) and velocity (Vsh) of the 
inner shell are found to be strongly correlated through conservation 
of energy and momentum. Processes such as turbulent mixing and 
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reverse shocks cause the thickness (Ai?sh) of the inner shell to re- 
main constant once it has begun to coast (defined to begin at radius, 
/?sh,o)- A detailed description of the formation of the inner shell as 
well as emission methodology used in this work can be found in 
paper 1. The velocity of the inner shell that we found as a fit pa- 
rameter in paper I was approximately 10% of our theoretically pre- 
dicted value. Upon revisiting our derivation of R band luminosity 
used in paper I, we found that there was a missing factor of tt and 
with this small correction the inner shell velocity used in this work 
is now in full agreement with the theoretical value (see eqn. 7 of 
paper I). 

The inner shell is fully parameterized by: temperature, mass, 
velocity, thickness and coasting radius. The mass and velocity are 
determined through conservation of energy and momentum, the 
derivation of the formulae used to describe these parameters can 
be found in paper I (equations 6 and 7 respectively). For our anal- 
ysis of each SLSN target considered in this work we use the same 
coasting radius for the inner shell, namely the best fit value of 
^sh,o = 4 X 10^"^ cm that was found in paper I. The only param- 
eter pertaining solely to the inner shell that we allow to be adjusted 
is shell thickness. 

A detailed look at the formation of the inner shell using a full 
hydrodynamic treatment would help to understand how physical 
parameters of the inner shell are aff'ected by changing initial condi- 
tions such as the time delay between SN and QN. 

4.2 The Envelope 

Fiducial values, found in paper I, for the shock speed (Vshock = 6000 
km s~^) and the outer edge velocity of the homologously expanding 
(v cx r) envelope (vsN,max = 4100 km s~^) are used for the work 
presented here. By fixing Vshock and vsN,max the variation in the time 
at which the cooling phase begins is uniquely determined by fdeiay 
Following the same methodology as paper I, a simple temperature 
profile is considered for the envelope where T(r) oc with j3 - 
0.2. In this work the progenitor star mass (M^) is a free parameter 
and the mass of the inner shell (M^h) is prescribed by conservation 
laws, thus the mass of the envelope is found simply as Mgnv = - 



4.3 Cooling 

Since the inner shell and the envelope are initially both parts of the 
same physical structure, namely the re- shocked SN envelope, they 
both start with same initial temperature (Tq). As the envelope and 
the inner shell expand they cool adiabatically, due to the diff'erence 
in geometries (the envelope expanding spherically and the inner 
shell with constant thickness) they follow diff'erent cooling profiles. 
The temperature of the envelope evolves as renv(0 ^ Tq and for 
the inner shell oc Tq r"^!^ (see paper I for derivation). 

As continuum radiation is emitted by the envelope an equiva- 
lent amount of thermal energy is removed, thereby conserving en- 
ergy. Since the cooling phase begins when the envelope is still rela- 
tively optically thick, radiative cooling starts at the outer edge of the 
envelope. Due to radiative cooling, over time a hot-cold interface 
progresses inward through the envelope. In our model once an outer 
layer of the envelope is radiatively cooled it no longer emits radi- 
ation. At the late stages of the dsQN light curve, radiative cooling 
will have caused the hot-cold interface to progress such that most 
of the envelope has been cooled. By this time the envelope has also 
expanded causing it to become more optically thin and thus mak- 
ing our optically thick radiative cooling approximation less valid. 



Table 2. List of SLSNe targets. 


Name 


Type 


Peak Magnitude 


Proposed Models 


SN 2005ap 


H-poor SLSN^ 


-22.7'' 


CSM^ PISN^dsQN'^ 


SN 2006gy 


Iln^ 


-22^ 


PISN^ CSM^ dsQN'^' 


SN 2006tf 


Iln^^ 


—20.8^ 


CSM*^ 


SN 2007bi 


Ic^ 


-21.3i 


CCj,PISN' 


SN 2008es 


II-L^ 


-22.2^ 


CSM^ 


SN 2008fz 


Ilni 


-22.3' 


none' 


PTF09cnd 


H-poor SLSN^ 


-22^ 


CSM^ PISN^ 


PTFlOcwr 


H-poor SLSN^ 


-21.2"^ 


CSM^ PISN^ 



^ Quimby et al. (201 1), Quimby et al. (2007b), ^ Leahy & Ouyed (2008), 
^ Quimby et al. (2005), ^ Smith & McCray (2007), ^ Woosley et al. (2007), 
g Ouyed et al. (2012), ^ Smith et al. (2008), ' Gal-Yam et al. (2009), J Young 
et al. (2010), ^ Gezari et al. (2009), ' Drake et al. (2010), Pastorello et al. 
(2010) 



In reality as the envelope becomes more diff'use radiative cooling 
will cause the whole envelope to cool rather than just the outer lay- 
ers and lead to an increase in the coohng rate of the envelope with 
time. The eff'ect of this can be seen in the envelope contribution to 
the dsQN light curve (plotted as a blue dotted line in Fig. 1): as the 
change in slope around day 170 is due to the envelope remaining 
hotter than physically likely. Since the envelope is optically thin 
when this coohng artefact becomes significant, in most cases lumi- 
nosity from the inner shell will already dominate the overall light 
curve rendering the problem moot. Rather than adding a more com- 
plicated radiative cooling law to our model we will simply note 
when the eff'ect is seen during our fits to observations. 



5 MODEL 

For this work we use the astrophysical modelling software shape, 
which allows us to construct the 3-D geometry of the dsQN sce- 
nario and perform radiative transfer (Steff'en et al. 2010). The tem- 
perature as well as the dimensions of the geometry are governed by 
the physics described in section 4. 



5.1 Radiative Transfer Parameters 

For this analysis the radiative transfer calculation follows the same 
methodology used in paper I in which for each frequency (v) an 
emission and absorption coefficient are specified. The emission co- 
efficient used for the envelope has the form 

;V = (1) 

where A is a multiplicative factor which is related to the underlying 
radiative process. For this analysis we chose to fix A over the filter 
passband of each studied SLSN. For the R band A = 5x 10^, V band 
A = lxlO^BbandA = 7xlO^,ubandA = 5 x 10^ The variations 
required in A allude to an emission mechanism more complex than 
our approximation described by eqn. 1. Assuming local thermody- 
namic equilibrium the absorption coefficient corresponding to our 
jy is found by means of the Planck function {By) to be 



(2) 



We as well include a Thompson scattering term to equation 2 of the 
form 



(3) 
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where ctth is the Thompson scattering cross-section and B is a mul- 
tiplicative factor which represents the fraction of scattered light 
that is not scattered back into the beam. For this work we set 
B = 5 X 10~^, which is the same as that used in paper I. Further 
details on shape and the radiative transfer calculation can be found 
in paper I. 



5.2 Fiducial Model Characteristics 

Before fitting the observed SLSNe light curves with the dsQN 
model we felt it would be informative to first explore the light 
curves of dsQNe in a more general sense. To this end we built 
a generic dsQN model with the following physical parameters: 

= 25Mo, To = 2.5 x 10^ K, A/?sh = 3 x lO^^cm and ^deiay = 15 
days. The light curve associated with this model is plotted as a 
red solid line in Fig. 1. The component of the light curve caused 
by emission from the envelope is plotted as a blue dotted line in 
Fig. 1 and the contribution from the inner shell is represented by 
the green dashed line. As can be seen in Fig. 1 the broad peak 
in the dsQN light curve is due to radiation emitted from the en- 
velope. As time progresses the envelope cools and becomes less 
dense causing a rapid drop in emission, this allows for radiation 
from the slowly coasting and thus slowly cooling inner shell to be- 
gin to shine through. The post-peak luminosity decay rate in the 
overall light curve of the dsQN is slowed due to emission from the 
inner shell. In the late stages of the dsQN hght curve there can exist 
a plateau due to emission from the inner shell. 

In this section we will investigate the eff'ect on light curve mor- 
phology of varying the physical parameters of a dsQN. For com- 
parative purposes in each of the panels in Fig. 2 the red sohd line 
denotes the same light curve as the overall hght curve seen in Fig. 1 
(red solid line). Then in each panel of Fig. 2 one physical parame- 
ter is adjusted to show its eff'ect on the light curve, the higher value 
is denoted by the blue dotted line and lower value by the green 
dashed line. Clock-wise starting from the upper-left panel of Fig. 
2; varying {20Mq, 25Mq and SOM©), varying Tq (2 x lO^K, 
2.5 X lO^K and 3 x lO^K), varying ^R,^ (2 x lO^^cm, 3 x lO^^cm 
and 4 x lO^^cm) and varying ^deiay (10 days, 15 days and 20 days). 

As seen in the upper-left panel of Fig. 2 the eff'ect of vary- 
ing is simply scaling the height of the peak, with higher mass 
yielding a higher peak. As the post-peak luminosity drops support 
by radiation from the underlying inner shell kicks in at the same 
time for each of the light curves. This convergence is caused by 
the fact that changing the mass of the envelope has no eff'ect on the 
emission from the inner shell. 

Changing Tq of the dsQN (as seen in the upper-right panel of 
Fig. 2) as well acts to scale the height of the peak, however the 
height of the peak is inversely proportional to Tq (due to the form 
of the emission coefficient, see eqn. 1). In the case of varying Tq 
the support from the inner shell on the overall light curve begins at 
diff'erent times because temperature directly eff'ects the amount of 
radiation that can be emitted by the inner shell. For higher Tq the 
luminosity of the inner shell is greater and its effect on the overall 
light curve can be seen sooner. 

As shown in Fig. 1 the eff'ect on the dsQN light curve of radi- 
ation from the inner shell can only be seen in the late stages once 
emission from the envelope has faded sufficiently. Thus as expected 
the lower-right panel of Fig. 2 shows that changing A/?sh only im- 
pacts the tail of the dsQN hght curve. All things being equal, a 
thinner shell implies higher density and thus more intense emis- 
sion from the inner shell. The eff'ect of increased density causing 



increased radiation can be seen in the lower-right panel on Fig. 2 
where the thinner A/?sh light curve has a more prominent plateau. 

The physical parameter that has the most dramatic eff'ect on 
light curve morphology is ^deiay» which can be seen in the lower-left 
panel of Fig. 2. A shorter ^^eiay implies that the SN envelope is more 
dense when it is re-shocked which translates to higher peak in the 
light curve. The shorter ^deiay also demands that the SN envelope is 
smaller when it is re-shocked, leading to significant adiabatic losses 
and a steep luminosity decay rate. The dsQN with a longer ^deiay can 
not achieve as high of a luminosity peak, however reduced adiabatic 
losses allow the light curve to remain at a higher relative brightness 
for a longer time when compared to a dsQN with a shorter ^deiay 
Another eff'ect of adjusting ^deiay is to change the time of peak lu- 
minosity relative to the time of SN explosion. A longer ^deiay dsQN 
will peak in brightness later in time due to the simple fact that the 
shock must traverse a larger envelope before it can break out. 

From Fig. 2 it is clear that changing the value of M^, Tq or 
^deiay ^11 eff'ect the height of the luminosity peak of the light curve 
of a dsQN. Our approximated emission coefficient (see eqn. 1) con- 
tains the free parameter A which as well directly eff'ects the amount 
of radiation emitted in our model at any given time. Thus we can 
use A to gain a better visualization of the eff'ect of changing each 
physical parameter (M^, Tq and ^deiay) on the shape of the light 
curve. In Fig. 3 we perform the same comparison as shown in Fig. 
2 with the addition that each light curve is rescaled using A such 
that the absolute magnitude of the light curve in the R band peaks 
at -22. By comparing the light curves in the upper two panels of 
Fig. 3 it is clear that neither changing nor Tq significantly af- 
fect the post-peak luminosity decay rate of the dsQN. As seen in 
the upper-left panel of Fig. 3 increasing the mass of the envelope 
marginally increases the breadth of the light curve peak but does 
not aff'ect the slope. Turning to the upper-right panel of Fig. 3 one 
can see that increasing Tq has the same eff'ect on the morphology of 
the light curve peak as increasing M^. The only diff'erence between 
adjusting and Tq is that the intensity of the radiation from the 
inner shell increases with Tq causing the plateau to vary directly 
with Tq. 

From Figures 2 and 3 we can see that varying Tq and ef- 
fectively only act as scaling factors, with minimal eff'ect on the lu- 
minosity decay rate of the light curve. Since the eff'ect is virtually 
the same, for all our fits to observed SLSN light curves we choose 
to fix Tq at 2.5 X lO^K and allow to vary. The lower-left panel of 
Fig. 3 displays the eff'ect of varying the delay between SN and QN. 
As was shown in Fig. 2 increasing ^deiay causes the peak luminosity 
to occur later with respect to the time of SN explosion. Another 
consequence of increasing ^deiay is to slow down the luminosity de- 
cay rate of the post-peak light curve of a dsQN. As we have found 
in this analysis the only parameter that we can adjust in order to 
change the luminosity decay rate in the dsQN model is ^deiay 



6 FITS TO OBSERVATIONS 

For this work we have fit the light curve of eight SLSNe using 
the dsQN model. The parameters that we adjusted to fit each set 
of observations were; fdeiay» A^sh and M^. Of these three free pa- 
rameters ^deiay and only significantly aff'ect radiation from the 
envelope while Ai?sh only aff'ects emission from the inner shell. As 
was shown in section 5.2, varying or Tq have the same eff'ect on 
the dsQN hght curve. Thus for simphcity we fixed Tq = 2.5 X 10^ K 
and allowed to vary between model fits. The best fit values of 
these parameters for each SLSNe surveyed can be found in Table 
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3. For each panel of Fig. 4 the observations are plotted with open 
circles and the best fit dsQN light curve is plotted as a red solid 
line. Fig. 5 displays each of the best fit dsQN light curves plotted 
on the same axes. From Fig. 5 it can be seen that for dsQNe in 
which an inner shell is formed (all but PTF09cnd and PTFlOcwr) 
the magnitude of the radiation from the inner shell is remarkably 
consistent. 



6.1 SN 2005ap 

The best fit progenitor mass for SN 2005ap was 28Mo. The high 
peak absolute magnitude and fast luminosity decay rate of SN 
2005ap required a short time delay (^deiay = 1 1.25 days) dsQN to fit 
the observations (see top row left panel of Fig. 4). Since the obser- 
vations of SN 2005 ap were only carried out for approximately 35 
days post-peak, the luminosity had not dropped enough to poten- 
tially see contribution from the inner shell. Thus we were unable to 
determine if an inner shell was present in the dsQN model of this 
SLSN. 



6.2 SN2006gy 

In the proof-of-principle paper (Ouyed et al. 2012) the fits to SN 
2006gy ignored a tt factor in the R band magnitude calculation, the 
fits presented here take this into account and thus are more accurate. 
Our fit to the observed light curve of SN 2006gy is displayed in the 
top row middle panel of Fig. 4. A SOM© progenitor star was used 
in our model and the broad peak of SN 200gy required a time delay 
of 17 days. The distinct plateau in the late stages of the fight curve 
was fit using a inner shell 2 x 10^^ cm thick. 

6.3 SN 2006tf 

Although the observations of SN 2006tf do not include the leading 
edge of the light curve, the slow luminosity decay rate demanded 
that the dsQN used in this fit have the longest time delay of all 
SLSNe modelled in this work (^deiay = 22.5 days). A progenitor star 
of 26Mq was used in our fit and the late stages of the light curve is 
aff'ected by contribution from an 4 x 10^^ cm thick inner shell. Our 
fit to the observed light curve of SN 2006tf is found in the top row 
right panel of Fig. 4. 

6.4 SN 2007bi 

A plot of our fit to the light curve of SN 2007bi is displayed in the 
middle row left panel of Fig. 4. A progenitor star of 25.5Mo was 
used in our model and the delay between SN and QN that best fit 
the shape of the SN 2007bi light curve was ^deiay = 18.5 days. The 
late stage plateau of the light curve was fit in our dsQN model by 
radiation predominantly from a2.5x 10^^ cm thick inner shell. 

6.5 SN 2008es 

The narrow peak of the light curve of SN 2008es was fit by dsQN 
of a 28Mo progenitor star with a time delay between SN and QN of 
12.5 days. At approximately 100 days after SN explosion the light 
curve of SN 2008es begins to be supported by luminosity originat- 
ing from a 3 X lO^^cm thick inner shell. A plot of our fit to the light 
curve of SN 2008es is displayed in the middle row middle panel of 
Fig. 4. 



6.6 SN 2008fz 

The combination of high peak luminosity and broad light curve 
made SN 2008fz require a 35Mo progenitor star and time delay 
of 20 days between SN and QN. Our dsQN fit to the light curve 
is found in the middle row middle panel of Fig. 4. The late stages 
of the SN 2008fz light curve is supported by luminosity from a 
5 X lO^^cm inner shell. 

6.7 PTF09cnd 

A plot of our fit to the light curve of PTF09cnd is displayed in 
the bottom row left panel of Fig. 4. The best fit parameters for 
PTF09cnd (^deiay = 17 days and - 31Mo) are similar to that 
of SN 2006gy with one exception. For PTF09cnd we found that the 
inner shell had to be nearly 10 times thicker than that of SN 2006gy 
(2 X 10^"^ cm, or approximately two fifths of the initial radius of the 
inner shell) in order to fit the late stages of the light curve. This 
diff'erence in inner shell geometry may be attributed to diff'erence 
in the abundance of the inner SN envelope and higher order mix- 
ing eff'ects related to specific dynamics of the QN/SN interaction. 
Another possible explanation for the reduced radiation from the in- 
ner shell at the late stages of the light curve of PTF09cnd is that 
the inner shell is cooling faster than the imposed T oc r^^l^ law. A 
larger cooling index would simply imply that the shell is expanding 
in thickness with time. A more detailed study of the hydrodynam- 
ics of the QN/SN interaction would help to further understand the 
formation process of the inner shell. The late stages of our dsQN 
model fit to PTF09cnd displays a change in slope which is an arte- 
fact of the radiative cooling approximation discussed in section 4.3. 

6.8 PTFlOcwr / SN 2010gx 

Displayed in the bottom row middle panel of Fig. 4 is the dsQN 
model fit to the observed light curve of PTFlOcwr. This SLSN dis- 
played the narrowest light curve of all those studied in this work 
and thus was fit with the shortest time delay (7.5 days) dsQN. A 
25Mo progenitor star was used in our best fit model and the fight 
curve of PTFlOcwr showed no sign of the presence of an inner 
shell. The late stages of our dsQN model fit to PTFlOcwr displays 
a change in slope which is an artefact of the radiative cooling ap- 
proximation discussed in section 4.3. 

7 DISCUSSION 
7.1 Spectra 

In the case of SN 2006gy there exists numerous observations that 
detail its spectral evolution (Smith et al. 2010). Unfortunately for 
other SLSNe, study of the spectra is limited by the amount of ob- 
servations available, which in many cases is one spectral obser- 
vation around the time of peak brightness. As discussed in paper 
I, the dsQN can display an Hof signature in which the underlying 
broad structure of the line is caused by thermally broadened emis- 
sion originating from the inner shell, and emission and absorption 
from the envelope result in a blue-side absorption feature. Fig. 6 
displays a schematic representation of the line of sight evolution of 
the Ha line in a typical dsQN. Unique to the dsQN scenario the 
wings of the Ha line are caused by thermal broadening rather than 
velocity broadening (see right-most Ha line in Fig. 6). With this de- 
scription the dsQN model can naturally account for the long-lasting 
breadth of the line that is inexpficable through other models (Smith 
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Table 3. Comparison of dsQN model parameters used to fit each SLSNe. 



SLSNe 


Time Delay (^deiay) 
[days] 


Mass ( M^) 
[Mo] 


Shell Thickness (Ai?sh) 
[10^3 cm] 


SN 2005ap 


11.25 


28 


N/A 


SN 2006gy 


17.0 


30 


2.0 


SN 2006tf 


22.5 


26 


4.0 


SN 2007bi 


18.5 


25.5 


2.5 


SN 2008es 


12.5 


28 


3.0 


SN 2008fz 


20.0 


35 


0.5 


PTF09cnd 


17.0 


31 


>20 


PTFlOcwr 


7.5 


25 


N/A 



et al. 2010). As the Hof line progresses through the hot component 
of the envelope a P Cygni profile is added to the structure of the 
line, this is displayed as the middle Ha line in Fig. 6. Finally the 
cold outer layer of the envelope adds absorption to the blue- side of 
the Ha line (seen in the left-most line plotted in Fig. 6). 

By classification SN 2005ap, PTF09cnd and PTFlOcwr do 
not display any prominent hydrogen lines. Although each of these 
SLSNe do show some indication of a weak Ha line, as mentioned 
in the supplemental material associated with Quimby et al. (2011). 
Unfortunately the poor signal to noise (S/N) ratio eliminates any 
possibility of studying these weak lines. The two spectral obser- 
vations of SN 2008fz each show a strong Ha line that resembles 
that of SN 2006gy (Drake et al. 2010). While the Ha lines display 
a broad underlying structure, the only spectral observation of SN 
2008fz that is at an epoch which would show significant absorption 
by the envelope has the red-wing of the Ha line cut oflF. Without the 
red- side of the line there is no way of comparing the observations 
to our predicted dsQN Ha line. The spectra of SN 2008es contains 
broad Balmer lines which display a P Cygni profile however the 
poor S/N ratio of the observations inhibits any study of these lines 
(Gezari et al. 2009). 

For our work in paper I we performed detailed analysis of the 
spectral evolution of the Ha line observed in SN 2006gy. The dsQN 
model was shown to provide a unique explanation for both the per- 
sistent broad structure of the line as well as the blue- side absorption 
feature seen in Ha of SN 2006gy. In this work we study the Ha line 
of two other SLSNe; SN 2006tf and SN 2007bi. For SN 2007bi the 
Ha line observed in the spectra 54 days after peak brightness (see 
upper left panel of Fig. 7) is mysterious as all other features of the 
spectra are consistent with a type Ic SN and thus should be free 
of hydrogen Gal- Yam et al. (2009). For comparative purposes we 
have plotted the Ha line observed in SN 2006gy (upper right panel 
of Fig. 7) from a similar epoch as the spectral observation of SN 
2007bi. The bottom panel of Fig. 7 is the dsQN model Ha line 
from the same epoch (from paper I). Although the S/N ratio of the 
SN 2007bi prevents the study of the delicate structure of the Ha 
line the large scale features are similar to that of a dsQN, namely 
the broad underlying structure and strong blue-side absorption. 

The spectrum of SN 2006tf has multiple observations at differ- 
ent epochs along its light curve. As noted by Smith et al. (2008) the 
Ha line observed in SN 2006tf resembles that of SN 2006gy as it 
contains a broad underlying structure that becomes more prominent 
with time. The upper left panel of Fig. 8 displays an overplotting of 
four observations of the Ha line of SN 2006tf. For comparison the 
upper right panel of Fig. 8 displays an overplotting of four spec- 
tral observations of SN 2006gy from similar epochs. The bottom 
panel of Fig. 8 displays dsQN model Ha lines from roughly simi- 
lar epochs as the observations of SN 2006tf and SN 2006gy. In the 



dsQN model Ha line and that of SN 2006tf and SN 2006gy there 
exists a broad component to the line as well as a blue-side absorp- 
tion feature that diminishes in strength over time. 

The diflFerences in the Ha lines of SN 2006tf and SN 2007bi 
are consistent with the diflPerent ^deiay of the dsQN description of 
these SLSNe. In the case of SN 2007bi the shorter time delay 
(^deiay = 18.5 days) implies that the envelope is still dense at the 
time of spectral observation and thus there exists a strong blue-side 
absorption feature. As for SN 2006tf the time delay is significantly 
longer (^deiay = 22.5 days) therefore the envelope is much more dif- 
fuse during the spectral observations and the blue-side absorption 
feature is weaker. 

7.2 The Envelope 

The mass of the envelope as well as the initial shock temperature 
both act only to scale up or down the peak of the dsQN light curve. 
We found that the only parameter capable of significantly changing 
the luminosity decay rate is the time delay between SN and QN. 
All SLSNe with X-ray observations show remarkably quiet X-ray 
production. This is expected in the context of a dsQN due to the 
cold outer layer of the envelope which would act as an efficient 
absorber of high energy radiation. The formulation that we used for 
radiative transfer coefficients (see eqn. 1 and 2) provided a good 
fit to each of the SLSN light curves that we studied. This could 
be an indication that the continuum emission mechanism for these 
SLSNe resembles recombination which has a similar form as eqn. 
1 . Further study of SLSNe in which the light curve was observed 
in several passbands would help to determine the true temperature 
dependence for the emission coefficient and help in understanding 
the radiation mechanism. 

7.3 The Inner Shell 

Luminosity from the inner shell only begins to aflTect the shape of 
the light curves well past peak luminosity (typically once the lumi- 
nosity has dropped to 19.5 absolute magnitude). The slowing of 

the luminosity decay rate, or plateauing of the light curve in the late 
stage is caused by radiation from the inner shell shining through the 
diflfuse envelope. 



8 CONCLUSIONS AND PREDICTIONS 

We have shown that the dsQN scenario can be used to explain the 
light curves of all eight SLSN targets studied. In the context of 
dsQNe, progenitor stars ranging between 25-35 M© provide ample 
energy to power the large radiated energy budget of SLSNe. 



© 0000 RAS, MNRAS 000, 000-000 



8 M. Kostka et al. 



We found that the physical parameter with the greatest impact 
on dsQN hght curve morphology was the time delay between SN 
and QN. Shorter time delay dsQN yield a peak magnitude that is 
higher and a faster luminosity decay rate (narrower light curve). 
While for longer time delays the peak magnitude is lower and the 
light curve is broader, A variation in time delay in the dsQN de- 
scription provides an explanation for the wide variety of SLSN light 
curve morphology. From our analysis we found that for shorter time 
delay dsQN the inner shell may not be formed. The implication be- 
ing that the energy that would go into forming the inner shell may 
instead be lost to pressure- volume work, however further study of 
the dynamics of the QN-SN interaction must be undertaken. 

We have also examined the singular Ha spectral line pro- 
file found in three diff'erent SLSNe observations (SN 2006gy, SN 
2006tf and SN 2007bi). The broad structure of the line is accounted 
for by thermally broadened emission from the inner shell, while 
the intermediate peak and blue- side absorption feature are due to 
contribution from the envelope. We found that the evolution of the 
blue-side absorption feature in the Ha line of SN 2006gy and SN 
2006tf is consistent with diffiision of the envelope. 

Unique to the dsQN scenario is the fact that any core collapse 
SN that leaves behind a massive neutron star can in turn undergo a 
QN explosion. This is due to the fact that the conditions of the inte- 
rior of the progenitor star determine whether the neutron star could 
become susceptible to QN collapse. There is no correlation to the 
progenitor star envelope (for example whether or not hydrogen is 
present), thus we expect a wide variety of types of SN can become 
super-luminous due to re-brightening via a QN collision. 

A distinguishing feature predicted by the dsQN model is a 
unique chemical abundance caused by the spallation of the SN en- 
velope by the QN ejecta (Ouyed et al. 2011a). Recent observations 
have suggested plausible signs of the QN were found in Cas A 
(Hwang & Laming 2012). As described in Ouyed et al. (201 la) the 
layer of the SN envelope undergoes spallation by the QN ejecta de- 
pends on the density of the envelope at the moment of impact. For 
shorter time delay QN the Ni layer would be destroyed in favour 
of the production of sub-^^Ni elements such as; Ti, V, Cr and Mn 
(Ouyed et al. 2011a). This process should also lead to mixing that 
would cause Fe to be found in the outer regions of a dsQN remnant 
(Ouyed et al. 2011a). The formation of "^"^Ti at the expense of ^^Ni 
will play an important role in the very late stage (t ~ 1000 days) 
luminosity of dsQN. For longer time delay QN the spallation pro- 
cess occurs in the carbon & oxygen layer of the SN envelope which 
leads to a unique chemical signature of the dsQN. Spallation of the 
carbon & oxygen layer would lead to an over- abundance of lithium 
in the dsQN remnant (Ouyed 2012). 

A dsQN is expected to emit two bursts of X-rays. The first X- 
ray emission event would occur when the shock from the original 
SN breaks out of the stellar envelope and the second analogously 
occurs for the QN shock break-out of the SN envelope. If the time 
delay between SN and QN is short then the X-ray bursts could in 
fact be overlapped leading to a broadened X-ray light curve. How- 
ever if the time delay is long then there should be two distinct X-ray 
peaks. 



REFERENCES 

Amett W. D., 1982, ApJ, 253, 785 

Blinnikov S. I., 2008, in American Institute of Physics Conference Se- 
ries, Vol. 1016, Origin of Matter and Evolution of Galaxies, T. Suda, 
T. Nozawa, A. Ohnishi, K. Kato, M. Y. Fujimoto, T. Kajino, & 
S. Kubono, ed., pp. 241-248 



Bodmer A. R., 1971, Phys. Rev. D, 4, 1601 

Bucciantini N., Quataert E., Metzger B. D., Thompson T. A., Arons J., Del 

Zanna L., 2009, MNRAS, 396, 2038 
Chevalier R. A., 1982, ApJ, 258, 790 
Chevalier R. A., 2012, ApJ Lett., 752, L2 
ChevaHer R. A., Irwin C. M., 2011, ApJ Lett., 729, L6 
ChevaHer R. A., Irwin C. M., 2012, ApJ Lett., 747, L17 
Drake A. J. et al., 2009, ApJ, 696, 870 
Drake A. J. et al., 2010, ApJ Lett., 718, L127 
Fryer C. L., 1999, ApJ, 522, 413 
Gal- Yam A. et al., 2009, Nature, 462, 624 
Gezari S. et al., 2009, ApJ, 690, 1313 
Ginzburg S., Balberg S., 2012, ArXiv e-prints 

Grassberg E. K., Imshennik V. S., Nadyozhin D. K., 1971, Astro. & Space 
Sci., 10, 28 

Hwang U., Laming J. M., 2012, ApJ, 746, 130 

Itoh N., 1970, Progress of Theoretical Physics, 44, 291 

Kasen D., Bildsten L., 2010, ApJ, 717, 245 

Kasen D., Woosley S. E., Heger A., 2011, ApJ, 734, 102 

Kawabata K. S., Tanaka M., Maeda K., Hattori T., Nomoto K., Tominaga 

N., Yamanaka M., 2009, ApJ, 697, 747 
Keranen R, Ouyed R., Jaikumar R, 2005, ApJ, 618, 485 
Kiewe M. et al., 2012, ApJ, 744, 10 
Koning N., Ouyed R., Leahy D., 2012, ArXiv e-prints 
Law N. M. et al., 2009, PASP, 121, 1395 
Leahy D., Ouyed R., 2008, MNRAS, 387, 1193 
Miller A. A. et al., 2009, ApJ, 690, 1303 

Moriya T. J., Blinnikov S. I., Tominaga N., Yoshida N., Tanaka M., Maeda 

K., Nomoto K., 2012, ArXiv e-prints 
Niebergal B., Ouyed R., Jaikumar R, 2010, Phys. Rev. C, 82, 062801 
Ouyed R., 2012, ArXiv e-prints 
Ouyed R., Dey J., Dey M., 2002, A&A, 390, L39 

Ouyed R., Kostka M., Koning N., Leahy D. A., Steffen W., 2012, MNRAS, 
3043 

Ouyed R., Leahy D., 2009, ApJ, 696, 562 
Ouyed R., Leahy D., Niebergal B., 2007a, A&A, 473, 357 
Ouyed R., Leahy D., Niebergal B., 2007b, A&A, 475, 63 
Ouyed R., Leahy D., Ouyed A., Jaikumar R, 201 la, Phys. Rev Lett., 107, 
151103 

Ouyed R., Leahy D., Staff J., Niebergal B., 2009, Advances in Astronomy, 
2009 

Ouyed R., Staff J., Jaikumar R, 201 lb, ApJ, 729, 60 
Pan T, Loeb A., Kasen D., 2012, MNRAS, 2979 
Pastorello A. et al., 2010, ApJ Lett, 724, L16 

Quimby R., Castro F., Mondol R, Caldwell J., Terrazas E., 2007a, Central 

Bureau Electronic Telegrams, 793, 1 
Quimby R. M., Aldering G., Wheeler J. C, Hoflich R, Akerlof C. W., 

Rykoff E. S., 2007b, ApJ Lett., 668, L99 
Quimby R. M., Castro F., Gerardy C. L., Hoefiich P., Kannappan S. J., 

Mondol P., Sellers M., Wheeler J. C, 2005, in Bulletin of the American 

Astronomical Society, Vol. 37, American Astronomical Society Meeting 

Abstracts, p. 171.02 
Quimby R. M. et al., 201 1, Nature, 474, 487 
Rau A. et al., 2009, PASP, 121, 1334 

Smith N., Chomock R., Li W., Ganeshalingam M., Silverman J. M., Foley 

R. J., Filippenko A. V, Barth A. J., 2008, ApJ, 686, 467 
Smith N., Chornock R., Silverman J. M., Filippenko A. V, Foley R. J., 

2010, ApJ, 709, 856 
Smith N. et al., 2007, ApJ, 666, 1116 
Smith N., McCray R., 2007, ApJ Lett., 671, L17 
Staff J. E., Ouyed R., Jaikumar R, 2006, ApJ Lett., 645, L145 
Steffen W, Koning N., Wenger S., Morisset C, Magnor M., 2010, ArXiv 
e-prints 

Vogt C, Rapp R., Ouyed R., 2004, Nuclear Physics A, 735, 543 

Weidemann V, 1977, A&A, 59, 411 

Witten E., 1984, Phys. Rev D, 30, 272 

Woosley S. E., 2010, ApJ Lett., 719, L204 

Woosley S. E., Bhnnikov S., Heger A., 2007, Nature, 450, 390 



© 0000 RAS, MNRAS 000, 000-000 



Super-luminous Supernovae as Quark Novae 

Young D. R. et al., 2010, A&A, 512, A70 

Yungelson L. R., van den Heuvel E. R J., Vink J. S., Portegies Zwart S. K, 
de Koter A., 2008, A&A, 477, 223 



© 0000 RAS, MNRAS 000, 000-000 



10 M. Kostka et al 



-23 n 




I \ L \ \ \ \ \ \ \ \ \ \ \ \ \ \ 

50 100 150 200 

Time (days) 



Figure 1. Plotted as a red solid line is the R-band light curve of the dsQN. Radiation from the two components of the dsQN are also plotted. Radiation from 
the envelope is represented by the blue dashed line and that from the inner shell is denoted by the green dotted line. 
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Figure 2. In each panel the overlay of the R-band light curve of dsQN models with three different sets of physical parameters is plotted. For each panel one 
physical parameter is varied while the remaining parameters are held constant (values given in text). Top-left: Mass of envelope (M^). = 20Mq (green 
dashed line), = 25Mq (red solid line) and = SOM© (blue dotted line). Top-right: Initial shock temperature (Tq). Tq = 2x 10^ K (green dashed line), 
Tq = 2.5 X 10^ K (red solid line) and Tq = 3 x 10^ K (blue dotted line). Bottom-left: Time delay between SN and QN (tdeiay)- ^deiay = 10 days (green dashed 
line), fdeiay = 15 days (red solid line) and fdeiay = 20 days (blue dotted line). Bottom-right: Inner shell thickness (AR^h). A/?sh = 2 x lO^^m (green dashed line), 
ARsh = 3 X lO^^m (red solid line) and A/?sh = 4 x lO^^m (blue dotted line). 
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Figure 3. In each panel the overlay of the R-band light curve of dsQN models with three different sets of physical parameters is plotted. Each light curve is 
scaled using the radiative transfer parameter A such that the peaks reach the same absolute magnitude. For each panel one physical parameter is varied while 
the remaining parameters are held constant (values given in text). Top-left: Mass of envelope (Mi,). M^, = 20Mq (green dashed line), = 25Mq (red solid 
line) and = SOM© (blue dotted line). Top-right: Initial shock temperature (Tq). Tq = 2 x 10^ K (green dashed line), Tq = 2.5 x 10^ K (red solid line) and 
To = 3 X 10^ K (blue dotted line). Bottom-left: Time delay between SN and QN (tdeiay)- ^deiay = 10 days (green dashed line), fdeiay = 15 days (red solid line) 
and ^deiay = 20 days (blue dotted line). 
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Figure 4. For each panel the observed light curve is plotted with open black circles and the dsQN model is plotted as a solid red line. Time since the inferred 
SN explosion is plotted along the horizontal axis and the absolute magnitude in the observed band is plotted on the vertical axis. The best fit parameters used 
to generate these fits can be found in table 3. From left-to-right and top-to-bottom the panels represent: SN 2005ap, SN 2006gy, SN 2006tf, SN 2007bi, SN 
2008es, SN 2008fz, PTF09cnd and PTFlOcwr. 
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Figure 5. The best fit dsQN model light curve for each of the SLSNe studied in this work are plotted on the same axis. 
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Figure 6. The three Hof emission lines plotted here represent the Ha line at three different stages of its line of sight evolution. In the background a not-to-scale 
cartoon representation of the physical structure that contributes to the given emission line. The far right Hof line is thermally broadened emission from only 
the inner shell. The broad structure of the middle Hof line is emission from only the inner shell, while the hot envelope contributes a P Cygni profile on top of 
the broad structure. The far left Hof line is emission from the entire dsQN along the line of sight. The radiatively cooled outer layer adds increased blue-side 
absorption to the Ha line. 
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Figure 7. Comparison ofUa line from spectra of SN 2006gy, SN 2007bi and the dsQN model. Top-Left: Plotted is the observed Ua hne of SN 2007bi approx- 
imately 54 days after peak luminosity, data from Gal- Yam et al. 2009. Top-Right: The Ua line from the spectrum of SN 2006gy observed at approximately 50 
days post-peak luminosity, data from Smith et al. 2010. Bottom: The dsQN model Ha line 50 days post-peak luminosity, from paper I. 
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Figure 8. Comparison of evolution of Ha line from SN 2006tf spectra, SN 2006gy and our dsQN model Ho- line evolution. Top-Left: Plotted is an overlay of 
spectral observations of the Ho- line observed in SN 2006tf (data from Smith et al. 2008). The observations are from; 32 (black solid line), 41 (blue dash line), 
66 (red dotted line) and 95 (green dash-dot line) days after the first observation. Top-Right An overlay of Ha spectral lines from SN 2006gy is plotted, data 
from Smith et al. 2010. The Ho- lines plotted are from days; 96 (black solid line), 125 (blue dashed line), 154 (red dotted line) and 179 (green dash-dot line) 
after the inferred SN explosion date. Bottom: Plotted is an overlay of the dsQN model Ha line that was used to fit the observations of SN 2006gy in paper 
I. The Ho- lines plotted are from days; 96 (black solid line), 125 (blue dashed line), 154 (red dotted line) and 179 (g''^^(g)'J||5\]'|j{^^5^^^^§'g(j9fs 



